First-order formalism for scalar field in cosmology 
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We present a general procedure to solve the equations of motion for cosmological models driven 
by real scalar fields with first-order differential equations. The method seems to have great power, 
since it works for closed, flat or open space-time, for scalar fields with both standard and tachyonic 
dynamics. We illustrate the procedure solving several examples which model situations of current 
interest to modern cosmology. 
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| Dark energy is some necessary but yet very mysterious type of energy. It is introduced to respond for about 3/4 
of the physical contents required to fuel the accelerated expansion of the Universe Q - see also Q for more recent 
investigations. The unusual properties of dark energy defy the interested community, which has introduced a diversity 
of mechanisms and models to explain its main characteristics - see, e.g., the recent works 0- A lesson to learn from 
the many investigations is that the possibility to describe properties of the dark energy in cosmological scenarios 
is related to the presence of scalar fields, which are simple, dynamical and allow for several distinct mechanisms of 
00 ■ current interest to modern cosmology. We can also include an exotic fluid which behaves as the Chapligin gas Q. 
£NJ | In this case, the fluid has a bizarre equation of state, making the pressure to depend on the negative of the inverse 
■ energy density, p ~ — 1/p. The Chapligin fluid has interesting non-standard features some of them bridging non- 
relativistic fluid dynamics behavior into brane-like properties of relativistic Nambu-Goto actions. In cosmology, an 
interesting property of the Chapligin fluid is that it may interpolate between different phases of the cosmic evolution, 
from dust-dominated to de Sitter universe, passing through intermediate phase which mix cosmological constant and 
stiff matter, where p = p, and this has inspired many recent investigations. 

An important difficulty related to dark energy is intrinsic to cosmology, which is governed by Einstein's equation. 
Even under the presence of very important assumptions, which lead to the Friedmann-Robertson- Walker (FRW) 
model, one still gets Einstein's equation in the form of two non-linear differential equations which are in general very 
hard to solve. In the presence of scalar fields, other equations should be added, the equations of motion for the several 
degrees of freedom introduced by the scalar fields present in the model. For this reason, investigations concerning the 
identification of general properties of the equations of motion in the presence of scalar fields have gained increasing 
, importance, as one sees for instance in the recent investigations 0,0- 

In the investigations presented in [||, some of us have shown how to trade the equations of motion to first-order 
differential equations, for cosmological models of the FRW type, driven by scalar fields with standard or tachyonic 
dynamics, in closed, flat or open space-time. Some of these results can be seen as generalizations of the Hamilton- 
Jacobi formulation of inflation which appeared before in Ref. Evidently, the presence of first-order differential 
equations is important to inject new motivation into the subject. 

The procedure introduced in |(j is inspired in the fact that in FRW cosmology driven by a real scalar field <j>, both 
the Hubble parameter and the scalar field depend on time explicitly; however, if we write t = t(<f>), we can turn the 
Hubble parameter a function of the scalar field, H = H((f>). In Ref. [6| we have implemented this possibility with the 
inclusion of another function, W — W ((/)), in the form H = W ((/)). In flat spacetime, this possibility works very nicely, 
if the scalar field potential can be written as a function of W in a specific manner, demanded by the equations of 
motion of the model under investigation. The key point here is that H = W((f>) can be seen as a first-order differential 
equation for the scale factor a = a(t) which defines the Hubble parameter H = a/ a, with dot standing for the time 
derivative. The use of H = W and the equation of motion for the acceleration gives another first-order equation, 
4> = —W$, with Wff, standing for dW/d(f>. Now, these two first-order equations and the Friedmann equation demand 
that the potential V = V(<fi) has the form V = 3W 2 /2 — This procedure is similar to the Hamilton- Jacobi 

formulationpresented in Chapter 2 of Ref. p| . 

In Ref. |3 we have been able to go beyond flat space-time, writing first-order differential equations for all the 
three distinct cases of spherical, flat and hyperbolic geometry. These results were later used to extend the first-order 
formalism to the braneworld scenario of the Randall- Sundrum type 0, which is driven by warped geometry with 
a single extra dimension of infinity extent |lfj. The braneworld side of the problem has led us to establish direct 
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connection with former investigations, both in flat 0] an d curved space-time 0,0] - see also Refs. |?l ITil ITE1, Il6j . 
which introduce several considerations concerning supersymmctry and the first-order formalism for brane and for 
cosmology. These findings are now used to further explore the first-order formalism put forward in Refs. |UE3- Below 
we consider models of FRW cosmology, driven by real scalar field which evolves with standard or tachyonic dynamics. 
In the case of standard evolution, we focus on the generality of the procedure and we illustrate the method with some 
examples of current interest to cosmology. For tachyonic dynamics, we also make specific progress, since we extend 
the result of 6], which was obtained in flat space-time, to the case of spherical or hyperbolic geometry. 



2. Standard dynamics 



To make the above reasonings sound, wc investigate the general problem which is described by the action 

d A x J—g ( -I R + £(</>, d^)) (1) 
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where R and <f> stand for the scalar curvature and scalar field, respectively, and we are using AttG — 1. As usual, we 
use g to represent the determinant of the metric tensor g^ , which can be obtained by the line element 

ds 2 = dt 2 - a 2 (t) ( + r 2 dQ 2 ) (2) 
\ 1 — kr z J 

In this expression, the last term describes the angular portion of the three-dimensional space, a(t) is the scale factor, 
and k is constant: k = 1,0, or — 1, for spherical, flat, or hyperbolic geometry, respectively. 

The energy- momentum tensor has the form T^ v = (p, —p, —p, —p), where p and p represent energy density and 
pressure, respectively, which can be expressed in terms of the scalar field described by the Lagrange density C. The 
above action leads to Einstein's equation, which allows obtaining 

- = -\(p + ty (3b) 
The equation of motion for the scalar field depends on C{4>, d^cj)); for standard dynamics we have 



c = -d^d^~v{4>) (4) 



and the equation of motion for <j) has the form 



3Hj> + %■ = (5) 



The energy density and pressure which account for standard dynamics are given by 



P = \tf + V, P=\fi-V (6) 



We use Eqs. © and J3J to write 
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H 2 = -tf+o V -— ( 7a ) 

(7b) 

The set of Eqs. (JjjJ and Q constitutes the equations we have to deal with to solve the problem specified by the 
potential V ((/)). 

To trade this set of Eqs. J5J and J7J) to first-order differential equations, we follow Refs. HG3. The first step is to 
make the choice 



H = W + akZ 



(8) 
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which leads to the equation 

4> = -W4, + f3kZ^ (9) 

where W = W((f>) and Z — Z(4>) are in principle arbitrary functions of the scalar field, and a and f3 are real parameters. 
These equations are first-order differential equations, and they solve the set of equations (JSJ and for the potential 

V = | (W + akZf + ~ (fikZ^ - Wt)[(3a + 2/?)fcZ + W+] (10) 

and the constraint 

W^Z^ + - 2akZZ^ - 20kZfZM - 2WZ^ = (11) 

We can make the choice a = 0. It gives 

H = W (12) 

and 

j> = -W4 + !3kZ^ (13) 

In this case we get the potential 

V=^W 2 + ^(fikZ^-W^)(20kZ^ + W 4> ) (14) 

and constraint 

+ Vf*Z H - WkZ^Z^ - 2WZ^ = (15) 

This choice leads us back to the investigations done in Ref. j^|, if we change Z^ to Z, as we have used in 6]. Another 
possibility is given by = 0, which leads to 

H = W + akZ (16) 

and 

4> = (17) 

In this case we get the potential 

V = ^{W + akZ) 2 -^W^akZ^ + W^ (18) 

and constraint 

W^Zj, + W+Zm - 2akZZ^ - 2WZ <j) = (19) 

To investigate the cosmic acceleration, instead of using the deceleration parameter q = —da/a 2 we introduce q = —q 
as the acceleration parameter, which has the form 

aa - H . . 



We use Eqs. © and © to get 



{W + akZf 1 ' 



For a = it gives 

W 



t | wMkZj -Wj) (22) 
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and for (3 = it becomes 



9=1 



W^Wi + akZt) 
(W + akZ) 2 



(23) 



and we see that there are several different possibilities to treat the acceleration. 

The choice a — or (3 — corresponds to different gauge choice - see, for instance, Refs. The case 

(3 = will be further investigated in the present work. It is described by the first-order equations l|16(l and (|17|l . and 
the scalar field model has the potential given by equation (|18|l . where the functions W{4>) and Z(<fi) are constrained 
to obey l|19|l . and the acceleration parameter has the form (|23[l . In this case, we notice that if Z = Z((f) is given, we 

can introduce W = W + akZ in order to rewrite (| 1 f>p and 117|l in the form H = W and 4> = —W<f, + akZ§. This pair 
of first-order equations has the very same form of the former pair 1|12|) and (|13[l . Thus, choosing the pair l|12l) and 
(|13[l or 116|) and (|17|l is just a matter of convenience. Evidently, we can also choose the pair JSJ and ©, but in this 
case one needs to set a + (3 ^ 0. 

We illustrate the subject with some examples. In the case (3 = 0, let us first consider the choice Z = <f>. In this case 
the constraint leads to 



- 2W - 2akcj) = 

We can solve this constraint with W = sinh(A</>) — akcj), for A = ±\/2. The scalar field has the form 

B 



A arctanh 



A 



( AB 



(24) 



(25) 



with B = V a 2 k 2 — 2, which requires that a 2 k 2 — 2 > 0. The potential is given by 



V = -sech 2 (Acj)) ~ [Acosh(A(j>) - ak] 



— cosh(A0) + ak 



(26) 



and now the Hubble parameter gets to the form 



H 



2B(A + ak) tanh(^t) 
(A + ak) 2 -B 2 tanh 2 (^) 



(27) 



The cosmic acceleration is given by 



q = 



A [(A + ak) 2 + 2 tanh 2 (4^)] sech 2 (4^) 
+ 4 (A + ak) tanh 2 (4fi) 



(28) 



It is plotted in Fig. 1 for for some specific values of a and k. We see that for k — — 1 the acceleration changes sign, 
nicely indicating that the cosmic evolution may evolves from deceleration to acceleration. 



k= 1 
k = -l 



FIG. 1: Plots of the acceleration II28H for a — 2, and k = ±1. 
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We take the same model, but we now choose Z — W. The constraint equation reduces to 

- (1 + ak)W = 



(29) 



We can solve this constraint with W = Acosh(Bcj)), com B = \/T+ ak. The scalar field is given by 



, ( AB2 
tanh £ 

V 2 



and the potential has the form 



The Hubble parameter is 



V = ^AB 2 cosh 2 (B4>) - ~(1 + 3ak)AB 2 sinh 2 (£<?!>) 



and the cosmic acceleration gets to the form 



AB 2 



tanh" 



tanh (^fit) 



4 tanh' 



(^*) 



[tanh 2 (^t) + l]" 

which is plotted in Fig. 2, showing a new scenario, well distinct from the former one. 



(30) 



(31) 



(32) 



(33) 




FIG. 2: Plots of the cosmic acceleration for A = ak = 1. 



3. Tachyonic dynamics 

We now turn attention to the case of tachyonic dynamics. The importance of tachyon in cosmology is inspired in 
String Theory, and can be seen by the recent works 0] . In this case, the Lagrange density is given by 



The scalar field now evolves differently, according to 

i+(l-j> 2 )(3Hj>+¥f) =0 

and the energy density p t and pressure pt are changed to 



Pt = V 



1 



Pt = -V^l 



(34) 



(35) 



(36) 
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We use Friedmann's equation and the equation for the acceleration © to write, in the presence of curvature 

H 2 = \^=~- 2 (37a) 
i 2 k 

H = 7==V + -= (37b) 

a 

We extend the first-order formalism of Ref. |(| to tachyonic dynamics in the presence of curvature as follows. We 
consider 



This choice allows writing the equation 



In this case, the potential has the form 



H = VW 2 - akZ 2 (38) 



(39) 



3Wy/(W 2 -akZ 2 ) 



V = - W 2 



■ n 1/2 

i _ 4_ W 2 /W 2 



(40) 



9 W 2 - akZ 2 
and the constraint becomes 

3W 3 Z - SakWZ 3 - 2W 4> Z ct> = (41) 

The acceleration parameter is given by 

2 WW* - akZZ^ 
h 3 W (W 2 -akZ 2 ) 2 { ' 

In the slow-roll approximation [l^, one uses that </> and 2 are very small. Here it makes no sense to neglect <j>, but 
if we consider cp 2 very small, we can rewrite the potential H40|) in the much simpler form 

3 1 Wl 

V = ± W 2-L V l± (43) 

2 3W 2 -akZ 2 1 ' 

and this is the result we could obtain if one uses (|38[) . for scalar field evolving in accord to the standard dynamics. 
This result is very natural: in the tachyonic Lagrange density l|34|) . if we expand the square root up to first order in 
d^d^cf), we get to standard dynamics. 

To illustrate the general situation, let us consider some examples. We firstly notice that the case Z = 1 is trivial, 
and so we consider Z = W. This case leads us to the constraint 

3(1 - ak)W 4 - 2W<p 2 = (44) 

which is solved by 

W = ± , =- (45) 

y/Z(l-ak)<j> 



In this case we get 



and the potential 
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The Hubble parameter is H = 1/i, and this leads to q = 0. 
Another possibility is Z = W 2 . Here we get the constraint 



which is solved by 



This gives 



and 



The Hubble parameter is 



and now we get 



3(1 - akW 2 )W 4 - AW^ 2 = 



W = ± 



^Aak + 3</> 2 
V3, 



V 



H 



2V6 



Aak + 30 2 



21 



Aak + 1 2 



Aak + 1 2 



2t 2 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



We plot this acceleration parameter q in Fig. 3. The new scenarios show that for tachyons the acceleration goes 
asymptotically to 1/2, but the cosmic evolution is very similar to the case shown in Fig. 1. 



— k= 1 

— k = -l 



FIG. 3: Plots of the cosmic acceleration 153H for a = 1, and k — ±1. 



4. Summary 

We summarize this work recalling that we have extended the first-order formalism introduced in Ref. to describe 
FRW models, driven by a single real scalar field with either standard or tachyonic dynamics for generic spherical, flat 
or hyperbolic spatial geometry. The importance of the procedure is related not only to the improvement of the precess 
of finding explicit solution, but also to the opening of another route, in which we can very fast and directly write 
the Hubble parameter H once Z — Z(cf>) and W — W((j>) are given. As we have shown, the present investigations 
are of direct interest to modern cosmology, since they seem to open several distinct possibilities of investigation. In 
particular, in the case of tachyonic dynamics, we have built the full formalism, including the cases of closed and open 
geometries, which were not given in the former work [||. The main results show that we can use standard or tachyonic 
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dynamics to find interesting cosmic evolution, including the case where the cosmic acceleration changes sign, evolving 
from deceleration to acceleration. 

The interest in the subject will certainly broaden with the extension of the method to the case of several fields, 
since in this case we could find models in which one field can be used to affect the behavior of other fields, unveiling 
the possibility to control a given phase and to link different phases of the cosmic evolution. Another possibility is 
related to the inclusion of the Chapligin £| and other fluid contents, such as dust, radiation and stiff matter. These 
issues are now under investigation, and we hope to report on them in the near future. 
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